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Effects of scale coupling between mesoscopic slowly-varying envelopes of liquid-solid profile and 
the underlying microscopic crystalline structure are studied in the phase-field-crystal (PFC) model. 
Such scale coupling leads to nonadiabatic corrections to the PFC amplitude equations, the effect 
of which increases strongly with decreasing system temperature below the melting point. This 
nonadiabatic amplitude representation is further coarse-grained for the derivation of effective sharp- 
interface equations of motion in the limit of small but finite interface thickness. We identify a 
generalized form of the Gibbs- Thomson relation with the incorporation of coupling and pinning 
effects of the crystalline lattice structure. This generalized interface equation can be reduced to 
the form of a driven sine-Gordon equation with KPZ nonlinearity, and be combined with other 
two dynamic equations in the sharp interface limit obeying the conservation condition of atomic 
number density in a liquid-solid system. A sample application to the study of crystal layer growth 
is given, and the corresponding analytic solutions showing lattice pinning and depinning effects 
and two distinct modes of continuous vs. nucleated growth are presented. We also identify the 
universal scaling behaviors governing the properties of pinning strength, surface tension, interface 
kinetic coefficient, and activation energy of atomic layer growth, which accommodate all range of 
liquid-solid interface thickness and different material elastic modulus. 

PACS numbers: 81.10.Aj, 05.70.Ln, 68.55.A- 



I. INTRODUCTION 



Continuum theories have been playing a continuously 
important role in modeling and understanding a wide 
range of complex nonequilibrium phenomena during ma- 
terials growth and processing. For the typical example 
of liquid-solid front motion and interface growth, sharp- 
interface or Stefan-type models have been used in early 
studies to examine various solidification phenomena such 
as dendritic growth and directional solidification of either 
pure systems or eutectic alloys [l| . Recent focus has been 
put on the continuum phase-field approach, which has 
become a widely-adopted method in materials modeling 
not only due to its computational advantage as compared 
to atomistic techniques and also to the complex moving- 
boundary problems encountered in sharp-interface mod- 
els, but also due to its vast applicability for a wide vari- 
ety of material phenomena including solidification, phase 
transformation, alloy decomposition, nucleation, defects 
evolution, nanostructure formation, etc. [2tS|- 

These continuum methods are well formulated for the 
description of long wavelength behavior of a system. To 
incorporate properties related to smaller-scale crystalline 
details which can have significant impact, additional as- 
sumptions or modifications are required. Examples in- 
clude the consideration of lattice anisotropy for surface 
tension and kinetics Q, and the incorporation of sub- 
sidiary fields describing system elasticity [5, 6], plasticity 
Q, or local crystal orientation Q in phase-field mod- 
els. However, effects associated with the discreteness of a 
crystalline system, such as the atomistic feature of lattice 
growth, are usually absent due to the nature of contin- 
uum description. Efforts to partially remedy this in some 
previous studies include, e.g., the adding of a periodic po- 



tential mimicking effects of crystalline lattice along the 
growth direction in the continuum modeling of surface 
roughening transition 0-[H|: although both the form of 
lattice potential (usually assumed as a sinusoidal function 
as part of the sine-Gordon Hamiltonian) and the associ- 
ated parameters were introduced phenomenologically. 

More systematic approaches based on some fundamen- 
tal microscopic-level theories are needed for the construc- 
tion of continuum field models that incorporate crys- 
talline/atomistic features. One of the recent advances on 
this front is the development of phase-field-crystal (PFC) 
methods [l2l . [l3| , in which the structure and dynamics of 
a solid system are described by a continuum local atomic 
density field that is spatially periodic and of atomistic 
resolution; thus the small length scale of crystalline lat- 
tice structure is intrinsically built into the continuum 
field description with diffusive dynamic time scales. Both 
free energy functionals and dynamics of the PFC models 
can be derived from atomic-scale theory through classi- 
cal density functional theory of freezing (CDFT) and the 
corresponding dynamic theory (DDFT), for both single- 
component and alloy systems [13l - ll7| . Properties associ- 
ated with crystalline nature of the system, such as elas- 
ticity, plasticity, multiple grain orientations, crystal sym- 
metries and anisotropy, are then naturally included, with 
no additional phenomenological assumptions needed as 
compared to conventional continuum field theories. This 
advantage has been verified in a large variety of applica- 
tions of PFC, ranging from structural, compositional, to 
nanoscale phenomena for both solid materials II3~ 
HI] and soft matters [HHi]. 

An important feature of the PFC methodology is the 
multiple scale description it provides, as can be seen from 
its amplitude representation. The system dynamics is de- 
scribed by the behavior of "slow" -scale (mesoscopic) am- 
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plitudes/envelopes of the underlying crystalline lattice, 
as a result of the amplitude expansion of PFC density 
field in either pure liquid-solid systems [lij [26| or binary 
alloys [Til . [20l [27j . Note that in these amplitude equation 
studies, although most lattice effects have been incor- 
porated in the variation of complex amplitudes (mainly 
via their phase dynamics [HI), the spatial scales of the 
mesoscopic amplitudes and the microscopic lattice struc- 
ture are assumed to be separated (i.e., the assumption of 
"adiabatic" expansion). However, this assumption only 
holds in the region of slowly varying density profile ei- 
ther close to the bulk state or for diffuse interfaces, and 
hence is valid only at high enough system temperature. 
In low or intermediate temperature regime showing sharp 
liquid-solid or grain-grain interfaces, amplitude variation 
around the interface would be of order close to the lattice 
periodicity; thus the two scales of amplitudes vs. lattice 
can no longer be separated, resulting in the "nonadia- 
batic" effect due to their coupling and interaction. Such 
scale coupling leads to an important effect of lattice pin- 
ning that plays a pivotal role in material growth and 
evolution, as first discussed by Pomeau [28| and later 
demonstrated in t he p henomena of fluid convection and 
pattern formation [29l l30j . To our knowledge, these scale 
coupling effects have not been addressed explicitly in all 
previous phase-field and PFC studies of solidification and 
crystal growth. 

In this paper we aim to identify these coupling effects 
between mesoscale structural amplitudes and the under- 
lying microscopic spatial scale of crystalline structure, 
via deriving the nonadiabatic amplitude representation 
of the PFC model. What we study here is the simplest 
PFC system: two-dimensional (2D), single-component, 
and of hexagonal crystalline symmetry, as our main fo- 
cus is on examining the fundamental aspects of scale cou- 
pling and bridging that are missing in previous research, 
and also on further completing the multi-scale features of 
the PFC methodology. The explicit expression of the re- 
sulting pinning force during liquid-solid interface motion, 
and also its scaling behavior with respect to the interface 
thickness, are determined in this work, through the ap- 
plication of sharp/thin interface approach (given finite 
interface thickness) to the amplitude equations. This 
leads to a new set of interface equations of motion, in 
particular a generalized Gibbs- Thomson relation that in- 
corporates the pinning term and also its reduced form of 
a driven sine- Gordon equation. The pinning of the inter- 
face to the underlying crystalline lattice structure, and 
the associated nonactivated vs. nucleated growth modes, 
can be determined from analytic solutions of the interface 
equations for the case of planar layer growth. 



II. NONADIABATIC COUPLING IN 
AMPLITUDE EQUATIONS 

In the PFC model for single-component systems, the 
dynamics of a rescaled atomic number density field n(r,t) 



is described in a dimensionless form [12Ml4| 

dn/dt = V 2 [-en + (V 2 + ql) 2 n - gn 2 + n 3 ] + V • 77, 

(1) 

where e measures the temperature distance from the 
melting point, g — (S/B 31 ) 1 ' 2 /2 with B x proportional 
to the bulk modulus, and we have go = 1 after rescaling 
over a length scale R of lattice spacing. The noise field 
77 has zero mean and obeys the correlations 

{r} a (r,t)rf(r',t')) = 2T k B T6(r - r')6(t - t')5 aff (2) 

with a,/3 = x,y, z, where Tq is a rescaled constant de- 
pending on B x and R and T is the system temper- 
ature. 

To derive the corresponding 2D amplitude equations in 
the limit of small e, we need to first distinguish the "slow" 
spatial and temporal scales for the amplitudes/envelopes 
of the structural profile, i.e., 

X = e 1/2 x, Y = e 1/2 y, T = et, (3) 

from the "fast" scales (x,y,t) of the underlying hexago- 
nal crystalline structure. We then expand the PFC model 
equation (fT]) based on this scale separation and also 
on a hybrid ap pro ach combining the standard multiple- 
scale expansion [3 1[ |32| and the idea of "Quick- and-Dirty" 
renormalization group method [H (see Ref. QJ] for de- 
tails) . To incorporate the coupling between these "slow" 
and "fast" scales, which leads to nonadiabatic corrections 
to the amplitude equations, we use an approach based on 
that given in Refs. [29|, [30] which address front motion 
and locking in periodic pattern formation during fluid 
convection. 

Following the steps of standard multiple-scale analysis 
[3H [HJ , the atomic density field n can be expanded as 

3 

n = n (X,Y,T) + ^ Aj(X, Y, T)e iq ^' r + c.c, (4) 

where g° are the three basic wave vectors for 2D hexag- 
onal structure (i.e., the 3 "fast" -scale base modes) 

( \/3« y\ . / y\ 

Qi = -9o I — x + 2 1 > «2 = 9oy, <?3 = <7o I —x - - I , 

(5) 

and the slow scaled fields, including Aj (complex am- 
plitudes of mode q°) and no (real amplitude of the 
zero wavenumber neutral mode as a result of PFC 
conserved dynamics), are represented as power series 
of e: Aj = e^ 2 A\ 1/2) + eAf + e 3 / 2 ^ 3/2) + • • • = 
E~=i £" l / 2 4" /2) , n = eV2 n (V2) + m (D + e 3/2 n (3/2) + 

••• = E™=i em/2 ™o n/2) - Note that in E q- © ni s ner 

harmonic terms have been neglected. 

From Eqs. ([3]) and (f4]) as well as the substitutions 
dx( y ) ®x(y) + el ^ 2 dx(Y) an d d t — > edx, we obtain the 
following expansion for the PFC equation (JTJ) in the ab- 
sence of noise: 
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ST, 



[Cn + V 2 (gn 2 n 3 )]= ed T n - £ V 2 ^ + £ 



3=1 



+ c.c. 



j/fc=l 3=1 



(6) 



where £ = e\ + eV 2 — V 2 (V 2 +g 2 ) 2 is the linear operator 
in PFC, [■ • - ] s refers to the slow-scale expansion to all 
orders of e, and V s = (dx, dy), V 2 = d\ + dy, and C s 3 = 
eVg+e 1 ' 2 (2iq Q J • V s ) — q 2 are slow operators. In Eq. ([5]), 
J 7 ;, is the slow-scale correspondence of the effective free 
energy T given below [with (V 2 + 2iq® ■ V) replaced by 
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{C s , + ql) and (V 2 



T = / dr U- 



3=1 



- g 2 ) replaced by (eV 2 

3 

-3n 2 -2 5 n )El^f 

3 = 1 



?o 2 )]: 



•2i^-V)^| 



-(6n -2p) ( JjA,-+c.c. 

u =1 



[(v 2 



3=1 

+ 6Ei^i 2 i^ 

3<fe 
1 3 1 

o5 n o + 



. (7) 



which is the same as the previous amplitude expansion 
result Q H, HI]; also f° r other variables f Pjk (j,k — 
1, 2, 3) and f Pi (i = 0, 3) related to higher harmonics, 



AfA% 



f Pl = 3q 2 [(6no - 2< ? )A 1 A| + 3 
/ P2 = 4g 2 [(3n - .g) A 2 + A*] , 

/ P3 = 3(? 2 [(6n - 2g)A* 2 A 3 + 3 (AxA 2 + A\A* 2 



fpo — ^% 



m 



2 A, 



a;a;- 



2g)A 1 A* 3 + 3 (A\ 

f Pll =4q 2 [(3no - .9) A 2 + 6,M5A5] , 
/ P33 = 4g 2 [(3no - g)A 2 + ZA\A* 2 Az] , 
f P]k =2lq 2 A 2 Al (jjtk). (8) 

As in the hybrid method developed in Ref. [14], the 
amplitude equations governing Aj and uq can be de- 
rived from the integration of Eq. ([6]) over eigenmodes 
{e-<V},i.e, 



J^fLjy^^[Cn + W 2 (gn 2 



rx+X x dx^_ fV+ X v d]/_ \ r 
Jx X x Jy X y \ J -' 



V 2 ( 5 n 2 



e -iq° r v' = o 

3 

* 3 )L = 0, (9) 



where \ x = a = a /(\/3/2) and X y = \/3a = 2clq (with 
a the atomic lattice spacing for the hexagonal/triangular 
structure and ap = 2n/qo, as illustrated in Fig. [I}, which 
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FIG. 1. (Color online) Schematic of a hexagonal/triangular 
lattice, with atomic lattice spacing a = ao/(\/3/2), ao = 
2n/qo, and 3 basic wave vectors (j = 1,2,3) expressed 
in Eq. ©. 



are the atomic spatial periods along the x and y direc- 
tions respectively. Note that Eq. ([9]) can be also viewed 
as the combination of the solvability conditions obtained 
at all different orders of e in multiple-scale expansion [3] ■ 
In the limit of e — > 0, i.e., close to the melting tem- 
perature, the spatial variation of amplitudes Aj (X, Y, T) 
and no (X, Y 7 T) is of much larger scale compared to the 
atomic lattice variation scales and X y . Thus "slow" 
and "fast" length scales in the integral of Eq. © can 
be separated as in standard multiple-scale analysis (i.e., 
(X, Y) and (x,y) be treated as independent variables), 
leading to the Ginzburg-Landau-type amplitude equa- 
tions obtained in previous studies (la . l2o| : 

dA j/d t = C^^-q 2 ^ 0n /8t^ 2 ^-. (10) 

Note that to derive Eq. (fTO)) the long- wavelength approx- 
imation Cj — V 2 + 2iq® ■ V — g 2 ~ — q$ has been used as 
before. 

However, such assumption of scale separation would 
not hold when e is of larger value (still small but fi- 
nite, corresponding to low/moderate material tempera- 
ture far enough from the melting point). Although the 
amplitudes/envelopes still vary slowly in the bulk, the 
interface, either between liquid and solid states or be- 
tween different grains, could be thin or sharp, with its 
width comparable to "fast" lattice scales (e.g., of few 
lattice spacings). Thus functions of Aj and hq in Eq. 
(J6J) can no longer be decoupled from the atomic-scale os- 
cillatory terms e l t lq i~ mqk " l ' r (with I and m integers and 
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j,k = 1,2,3) in the integration of * dx' f y y dy' 
in Eq. ©• Using an approximation similar to that in 
Ref. [30l ] . we only keep the lowest-order coupling terms, 
i.e., terms coupled to lowest modes (with largest atomic 
lengths) along both x and y directions, including e lq i ' r , 



dAx/dt 



dA 2 /dt 



dA 3 /dt 



dno I dt 



where a projection procedure has been applied to address 
the noise term of the PFC Eqs. flTJ| and (0) Q, leading 
to zero mean of noise amplitudes rjj and tj , as well as 
the correlations (rji'Hj) = {voVj) — {VoVj) — an d 

(77,77*) = 2^T Q k B TS{r - r')5(t - t')S tJ , 
(rjgrfi) = 279 r fc B T ( 5(r - r')5(t - t')8 a ^ (15) 

with i,j = 1,2,3, a,j3 = x,y, and i?< = do = 
1/7 if assuming equal contribution from all eigenmodes 

{e lq j' v , 1}. In the above generalized amplitude equations 
(|li p -(|14 p . the integration terms explicitly yield the cou- 
pling between "slow" (for structural amplitudes or en- 
velopes) and "fast" (for atomic lattice variations) spatial 
scales, that is, the nonadiabatic corrections. The first 3 
coupling terms in each of Eqs. (|li p -(|13 p for the dynam- 
ics of complex amplitudes Aj are associated with lat- 
tice modes e %q i' r of length scale oq (= y/3a/2), while the 



corresponding to atomic layer spacing \/3a/2 — ao, and 
e*C'i — 9*«5 4 *" r ) with atomic layer spacing a/2 = ao/y/3. 
Couplings to higher modes, such as e 2lq i ' r (with length 
scale no/2), e^-^h^-r ( with length scale ao /V7), 
etc., are neglected. The following scale-coupled ampli- 
tude equations can then be derived from Eqs. © and 
©: 



I 

other 4 terms correspond to the coupling to e^ q i~ qk ^ j * k ' r 
modes with length scale ao/V3 (= a/2). The nonadia- 
batic effect for no dynamics is weaker, with only cou- 
plings to the 00/1/3 length scale given in Eq. (|T4|) . Note 
that in the bulk state of single crystal or homogeneous 
liquid, the amplitude functions f Pi , f Pjk ~ constants, and 
hence all the integrals in Eqs. (fl"T j) ~ (fT5)) are equal to zero; 
we can then recover the original amplitude equations (|10p 
without any nonadiabatic coupling, as expected. 



III. INTERFACE EQUATIONS OF MOTION 
WITH LATTICE PINNING 

To illustrate the important effects of scale coupling 
identified in above nonadiabatic amplitude equations, 
here we consider a system of coexisting liquid and 
solid phases, with the average interface normal direction 



ST 
SA* 
1 



AieA^ J x 

ST 



V+Xy 

dx' I dy 

>y 



— iqoV3x' ,f - 



iqo[^fx'+§y') »* iqo (^'-fl)') , <■* '90 (-^x'+f 2/') 



-4 



5A* 2 
1 



A^A^ Jx 



£C + A X 



y+x y 
dx' I dy' 
>y 



io (-^-x'+hy' 



if* AqoVSx' 1 f* —iqoVSx' . r ( -y^' + lf') j. iqo ( r x '+%Vj 

+ / Pll e + /p 33 e +./P2i e V ' +/P23 e V ' 



+ V2, 



=-% 



ST 

SA* 



An; Ay Jx 



X+Xx 



y+X y 
dx' I dy' 
'y 



f p -iq y' , f* Jlo(^x'+^y') , iqot^Lx'-^y') 
JP3 C <Jvn c "rJP33 L 



V3 „/ 3. 



Va, 



ST 
Snn 



\ x A 



x+K 



y Jx 



y+Xy 
d.r' I dy' 
'y 



f* oB iqoV3x' + f ^ e m{^x'+W) + fp Jqo(^x'-W) + , 



(11) 



(12) 



(13) 



V- 770,(14) 
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pointed along y. The extension of our derivation and 
results to other interface orientations is straightforward. 

In this case, the "slow" and "fast" scales parallel to 
the interface can be well separated, and to lowest-order 
approximation the amplitude equations (|ll[) -(|13 p are 
rewritten as 



dAj/dt : 



1 



A 



y+\ y 



v Jy 



-Q 2 o 



(v 2 + ^.v) 2 ^- + ^ 



V + >*y 



dno/dt = V : 



-7- / dy'f P] e^°y 



(16) 



V • 770 = V > + V • rjo 



(V 4 + 2,72V 2 ) n 



df_ 

dn 



V • r7 (17) 



( e -W : f or A t anc j a 3 ; e +i i°y' : for A 2 ), where \i = 
SJ^/Sn-o is a chemical potential of the system, and / is 
the bulk free energy density given from Eq. (|7|): 



3 3 
/ ={-e + 3n 2 - 2gn ) £ \A 3 \ 2 + | £ |^| 4 
i=i j'=i 



f 3 






+ c.c. j + 








1 3 , 1 4 




- 3^0 + l n o 



j<k 



(18) 



To derive the corresponding equations of motion for 
the interface with finite thickness £ (i.e., in the sharp/thin 
interface limit (3, 4]), we follow the general approach de- 
veloped by Elder et al. [|| with the use of projection 
operator method. In this approach, a small parameter e 
is introduced, which represents the role of the interface 
Peclet number, and the system is partitioned into two 
regions: an inner region around the interface, defined by 
— C < u < +£, and an outer region far from the interface 
(i.e., \u\ > £), where length £ scales as 1 -C £/£ -C e -1 , 
and u is the component of a local curvilinear coordinate 
in the interface normal direction. In this curvilinear coor- 
dinate (u, s), the two orthogonal unit vectors are defined 
as n = x sin 9 + y cos 9 (local normal of the interface) and 
t = dh/d9 (tangent to the interface), where 9 is the angle 
between n and the y axis; also 



V = nd u + i 



V 2 = d 



1 + UK 
K 



-d u 



if 



ud s n 

u ' l + un^ u ' (1 + un) 2 ~ (1 + uk) 3 

d9/ds. 



with the local curvature k 



A. Outer equations 

In the outer region which is far enough from the in- 
terface and close to the bulk states, the scale coupling 
term in Eq. (1161) can be neglected, and the slowly vary- 
ing amplitude fields Aj and no depend on rescaled spatial 
variables (eu/£,£s/£) and rescaled time e 2 t. Expanding 
the outer solution of amplitudes in powers of e, i.e., 



^out = j\0° 



f eA° ut + • • • , ng ut = 

..out ..out 1 ~~out 

/i — /x -r £fJ>i 



QOUt 



(20) 



substituting them into Eqs. ([To]) and ([TTJ), and using the 
rescaling given above, we find that at 0(1), 



df 



dA* 



and at 0(e), 



0, dn° ° Ut /dt = V 2 £L 



2,, out 



(21) 



df 



dA* 



= 0. 



dn° ut 



on 



V 2 ^, 



(22) 



out n ou t\ 

j . «o ) 
"|° ut " refers 



where "|g ut " refers to replacing (A,-, no) by (A° 
in the derivative df /dA* or df/dno, and 
to the corresponding results up to 1st order of A° ut and 
ng ut . Note that if assuming the system to be not far 
from a liquid-solid equilibrium state, Eq. (|2"Tj) of 0(1) 
yields the bulk equilibrium solutions of the uniform liquid 
(u = +00) or solid (it = —00) state 

Af u \u) ee A°(±oo), n°° ut («) = ng (±00), (23) 
with the corresponding equilibrium chemical potential 



B. Inner expansion and lattice coupling effect 

For the inner region (— £ < u < +£), the amplitudes 
and chemical potential can be also expanded as 



Af = A° + sAj 
/x in 



n™ = njj + en a + 
. 



(24) 



Due to the presence of interface at u = 0, the amplitudes 
are expected to vary rapidly along the normal direction 
n but slowly along the arclength s of the interface, lead- 
ing to the rescaling (U, S) = (u/£,es/£). Considering 
small interface fluctuations and noise amplitude, we as- 
sume that k — £«/£, 9 — e9, r]j — efjj, and rj = efj . 
Thus from Eq. (JTSJ we have V • rj = edufft/£ + 0(e 2 ), 
V 2 = (9 2 + ekdu)/^ + 0(e 2 ), and (V 2 + 2iq] ■ V) 2 = 
(V 2 + 2iq° j ■ V)g + e(V 2 + 2iq) ■ V)\ + 0(e 2 ), where for 
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J = 1,3, 

(v 2 + ^.v)^=(^_^ % ) 2 /e 4 , 

(V 2 + 2iq] ■ V)l =[Ti2V3q (sm6 duffi - iqotfu) 

+kd u (2d 2 u - 3q 2 e) T 2V3q 2 edud s 
-2t qo aa±V3ds)d 2 u \/^, (25) 

and for j = 2, 

(V 2 + 2zq 2 ° • V)" = (fl» + 2iq^du) 2 /t\ 

(V 2 + 2iql ■ V)\ =2Rdl{d u + 2iq d)/e- (26) 

To address the time relaxation of system in the inner 
region, as usual we use a coordinate frame co-moving 
with the interface at a normal velocity v n (s) = e^v(S) + 
0(e 2 ), and hence d t -» d t ~v- V = -evd v + 0(e 2 ). The 
inner expansion of the nonadiabatic amplitude equations 
(HU) and HZI) can then be given by: For 0(1), 



df 
i)A\ 



f) 2 



(V 2 + 2zqO.V)^o 

kd 2 (d 2 + 2q 2 e) nl + % 



0. 

= 0, 



(27) 



giving the equilibrium chemical potential /i cq = fio 
d 2 u {d 2 u + 2q 2 )nl + (df/dn )\ ; At 0(e), 



- g 2 (V 2 + 2z q °.V)^ J 



-% 2 E 



a 2 / 



fc=i 



a 2 / 



dA^dAl 



9 2 / 



n -9o 2 (V 2 + 2z^-V) 2 ^ 



Ju 



1 



1 



e Tiqo(U' + h) 



^^{di+qie)^ 

1 d 2 f 



4 d 2 {d 2 u + 2q 2 ^ 

d 2 f 



no 



-E 



dnodAj 



A: 



d 2 f 








(29) 



where we have assumed that the nonadiabatic scale cou- 
pling effects are of 0(e). For such nonadiabatic term 
appearing at the end of Eq. (|2"51) . X y — \ y /£, qo = qoS,, 
h = h/£ with h(x,t) the interface height, and we have 
used the transformation 



u = (y - h) cos 6 = (y - h)/y/l + {d x hf 



(30) 



with the lowest-order approximation u ~ y — h + 0(e 2 ) 
and dy ~ du + 0(e). 

Multiplying both sides of Eq. ([25} by 9c/ integrat- 
ing over J"_t dU (with £ = £/£ 3> 1), and then adding the 
results for all j = 1,2,3 and the corresponding complex 
conjugates, we obtain 



2{i E / _dU\duA° H \ 



-? 2 E 



dU 



(h dvAf - A*dun° 



a A K-p sin(q h + (p) 
d 2 f 



dA*8n 



dU^mduAf 



(31) 



where the boundary conditions 9jjA°(±£) = for any 
order k of derivative and ( 3> 1 have been used, 



/_^4e(i^°i 



es, \d V A 



-iqo£ 



{d 2 A\) (duAf) + (d 2 A°) (d v 



-4 



(32) 



-2(d 2 A°) (duAl*)}} 
with Si = Ss = 3<7q/2 and 62 — 0, and 

p e 1 ^-^ =2 j du [A\d u f v M%0 



+ATd u f P2 (A° j ,ri ) + A°d u f; 3 (A°,n ^3) 

with po > 0. Detailed derivation for this lattice coupling 
term can be found in Appendix [5J 

For Eq. (|29|) derived from the conserved dynamics of 
no, we need to adopt a Green's function method [4|. Sim- 
ilarly, two Green's functions are introduced, including 
G+(U,S;U',S') in the region < U < ( with surface 
S + closed at S = ±00, and G~{U, S\ U' , S') in the region 
— ( < U < with the corresponding surface 5_; they 
satisfy the equation 



a 2 / G ± (L7, 5; U', S') = S(U - U')S(S - S'), 



(34) 



with the boundary conditions G ± (t7 = U' = 0) = and 
d u G ± (U = ±C) = d w G ± (U' = ±C) = 0. Multiplying 
Eq. (|29| by G + (G _ ) and integrating over the corre- 
sponding region lead to 



dU' j> dS'v{S') G ± (U, S; U' , S') d uin ° 
o± Js± 

)\w=o± 

±£ ( ^ dU' <f dS'G ± d ul f]l l . (35) 
Jo ± Js± 



fn(U,S) ± <f dS' (G ± d w fi 1 - fnd w G^ 
Js+ 
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c.c. 



Further integrating Eq. (|35|) by J dU OuTIq and using the 
solutions of [see Eq. (|B1I) in Appendix [B] , we find 

vef_ dU [n§ - ng(±C)] 2 = - Ang#i (0, 5) - 
-E fdu(n d v Af -A*dv4) J£L 

+(dufa) ± { J _dU [n§-n§(±C)] 

-€ | C _ dtf K " «o(±0] [% u " %"(±0] , (36) 

where ng(d=C) = p°(C) for > and = rag(-C) for ?7 < 
0, 77o"(±C) - t?o(C) for C/ > and = fjjjf(-C) for U < 0, 
the miscibility gap 

An° = n°(C) - ng(-C) * n°(+oo) - ng(-oo) (37) 
due to £ 3> 1 in the inner region, and 



2 f c 



{dhnlY-qie {dunl 



(38) 



Also, the integration of Eq. (|29|) over (— £, £) yields the 
conservation condition for the inner solution 

- S^Ang = (duh)c ~ + £ [%"(0 " %"(-0] ■ 

(39) 

Combining Eqs. ff3T|) . ([36|) . and (|39|) and returning 
to the original unsealed coordinates (u, s), we obtain the 
following equation governing the normal velocity v n of 
the interface (given £ — > oo for the inner region) 



/+oo j 
du\2j2\duA j\ 
[ j 



g 2 [n°-n°(±oo)] ; 



= -<7 2 Ano^i(0, s) -aK-po sm(q h + tp) 

p OO 

-v n q 2 An° du[n° -n° (+oo)] + r] in , (40) 



where /xi = efi\ — p, m — /i cq + 0(e 2 ), the noise ryi n 
-q*f^ du[n° - ng(±oo)]^ - [J^du^mduAf 
c.c], and the surface tension cr is determined by a a 
ql& n , i.e., 

/-foo j 

+2% [(a 2 ^?) (a„<) + (<fr4g) (a u A° 
-2 (a 2 ^) (a u <)] + («) 2 - 9o 2 M) 2 }4i) 

Note that to derive Eq. ([^0]) , we have used the condition 



-j- OG 



<fo [ng - ng(±oo)] = 



(42) 



for a Gibbs surface to define the interface position u = 
[Hi- We find that this condition can also be derived at 
0(e), as shown in Appendix 151 



C. Results of interface equations 

To match the inner and outer solutions, we need to use 
the boundary conditions at u = ±£, i.e., 



Vi(u = ±(,s) = At ^ t (« = ±C,«), 



( ^rUic 



(43) 



and carry out the expansion of outer solution n° nt around 
the boundary. Based on the derivation given in Appendix 
IB1 from Eqs. (|4T)|) and (|B8I) we can obtain a modified 
form of the Gibbs-Thomson relation which incorporates 
the effect of coupling to the underlying lattice 



( v n = A - an -p Q sm(q Q h + tp) + r] v , (44) 



where 



(45) 



with 5/i = /i out - /i oq (= e/ii 11 '), and Ang = ng(+oo) - 
n[J(— oo) as defined in Eq. (1371) . which represents the mis- 
cibility gap given by the difference between bulk equilib- 
rium densities of coexisting liquid and solid states. Val- 
ues of Auq are small but always nonzero below the melt- 
ing point due to the first-order and metastability charac- 
ter of the liquid-solid transition. Also, Co is the kinetic 
coefficient determined by 



Co 



du\2Y,\duA°\ 



(±oo) 



(46) 

The noise term, rj v = rj; n + rj m (with rj m determined in 
Appendix [Bj , has zero mean and the correlation 



(r, v (s, t)r) v (s', t')) = 2DS(s - s')S(t - t'), 



(47) 



where D = i}q^T a k B T( , and ■& = di = t9 = 1/7 as in 
Eq. 05). 

Also the standard continuity condition for interface can 
be obtained from Eq. (j3"9")l and the matching conditions 
(see Appendix [B]), i.e., 



a o d5 V 
v n An = 



88 n 
du 



o+ 



(V5M) solid - (v^) Uq uid 



(48) 



where V^/x is evaluated at the location of moving solid- 
liquid interface. Finally to obtain the chemical poten- 
tial deviation at the interface from the outer solution 
8Aj = A° ut - A°° ut and Sn = n° ut - n°° u \ we need the 

which can be rewritten as 



_ A out — A 
lst-order outer equation 



df 



dA* 



where 



out 



2 df 



= 0, dSn /dt = V 2 — - 
on 



= V 2 5p, (49) 



" corresponds to the results of expansion up to 
1st order of 5Aj and Sng in the derivatives df /dA* and 
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df/dno. Note that from the equations (df/dA*)\° ut = 
(j = 1,2,3), each amplitude SAj can be expressed as a 
linear function of Sn , and hence Eq. (j49|) reduces to 
a diffusion equation of Suq with the effective diffusion 
constant depending on Aj° ut and 7ig° ut . 

The combination of Eqs. (|4"9]l. ([S]). and (j48]) yields a 
free-boundary problem, and can be reduced to the stan- 
dard form of sharp-interface equations if we neglect the 
lattice coupling term poshi(qoh + ip) in Eq. (|44[) . The 
incorporation of such scale coupling effect is analogous 
to the case of driven sine-Gordon equation describing the 
roughening properties of interface subjected to a periodic 
pinning potential [1(3, EH, or to the case of front lock- 
ing/pinning in fluid pattern formation [2^, Given 
v n = -du/dt ~ dh/dt/ll + idxh) 2 } 1 / 2 from Eq. O and 
k = V • h — — d%h/[l + (d x h) 2 ] 3 / 2 , for small local surface 
gradient d x h Eq. (|44|) can be approximated as 



Codh/dt = F + ad^h + 



\{d x hf 



Po sm(q Q h 



<f) + Vv 
(50) 

This has the same form as the (1+1)D version of the 
driven sine- Gordon equation introduced by Hwa, Kar- 
dar and Paczuski It is a variation of the sine- 

Gordon equation studied earlier by Nozieres and Gallet 
with an additional KPZ nonlinear term X(d x h) 2 /2 
Here Fq = A represents a thermodynamic driv- 
ing force determined by the chemical potential differ- 
ence Sfi at the interface y = h (i.e., u = 0), and terms 
ad 2 h — po sin(qoh + ip) can be derived from the sine- 
Gordon Hamiltonian. Compared to previous studies, our 
results given here in the PFC framework can determine 
detailed properties of the important parameters involved 
(including the kinetic coefficient £q, surface tension a, 
pinning strength po, and the driving force A), in partic- 
ular the explicit dependence on system temperature and 
elastic constants. Some example results will be given in 
the next section. However, it is important to note that 
while the above equation ([50)1 exhibits as a nonconserved 
form of interface dynamics, it is not complete and should 
be combined with Eqs. (|49| and (|48|) due to the condition 
of mass conservation required in a liquid-solid system. 



IV. APPLICATIONS TO CRYSTAL LAYER 
GROWTH AND PINNING 

To illustrate the important effects of nonadiabatic scale 
coupling on the dynamics of interface, we apply the in- 
terface equations of motion derived above to a simplified 
case of layer-by-layer crystal growth. The results, in par- 
ticular the different crystal growth modes of "continuous" 
vs. "activated" as well as the temperature and elastic- 
constant dependence of lattice pinning effect, can be used 
for examining the formation and evolution of more com- 
plicated surface/interface structures or patterns in fur- 
ther studies, the details of which will be presented else- 
where. For simplicity, in the following we consider the 
long wavelength limit of the average density field no and 



hence neglect the gradient terms of no in Eq. ([7]) for the 
free energy functional F (i.e., [(V 2 4- <? 2 )n ] 2 — > <Zo n o)j 
as such terms usually yield higher-order contributions 
to system properties [261 ] - The corresponding interface 
equations of motion given in Sec. IIII CI remain un- 
changed, although in Eq. (|4"Tj) for the expression of a 
the gradients of can then be neglected. 



A. Properties of interface parameters 

One of the most important parameters given in the 
above derivations is the strength of interface pinning 
force po. As determined by Eq. (f3"3")) . it depends on 
the details of liquid-solid equilibrium profiles A® and rig. 
These profiles are obtained by numerically solving the 
ID Oth-order amplitude equations given by (|27l) in an 
unsealed form: (V 2 + 2iq? • V)§A$ + (df/dA*)\ = 0, 
and d 2 (df/dno)\o — f° r the long wavelength limit of 
n®. We use a pseudospectral method in numerical cal- 
culations, and apply the periodic boundary condition by 
setting the initial configuration as 2 symmetric liquid- 
solid interfaces located at y = L y /A and 3L y /A. The ID 
system size L y perpendicular to the interface is chosen 
as L y — 2048Ay for all the results shown here, and a 
numerical grid spacing Ay = (27r/go)/8 is used. 

As shown in Fig. [2] (a), the pinning strength po in- 
creases with the decrease of system temperature (i.e., 
with the increasing value of e; see the inset), and also 
with the decrease of bulk elastic modulus B x for large 
enough po (> 10~ 14 ). This can be attributed to the phe- 
nomenon of sharper liquid-solid interface at lower tem- 
perature and smaller value of B x [see Fig. [5] (b)], since 
sharper interface leads to stronger scale coupling between 
microscopic crystalline structure and mesoscopic ampli- 
tudes, and hence larger pinning force; this is a funda- 
mental mechanism underlying the nonadiabatic deriva- 
tion given in Sec. [TXJ Thus one would expect that there 
might be a more universal relation between the pinning 
force and the interface thickness £, as can be derived from 
Eq. (|33[) governing po : Recalling that both amplitudes 
A® and rig are functions of scaled variable U = u/£ in 
the inner region (see Sec. IIII B[) . we rewrite Eq. (|33p as 



Po = 



due iqoU G(u/£) 



(51) 



where G = A\d u f^ +A° 2 *d u f P2 +A° 3 d u f^ Applying the 
residue theorem to the integral in Eq. (|51[) and assuming 
that within the poles (singularities) of G(U = u/£) in 
the upper-half complex plane, the one nearest to the real 
axis is given by U z — w z /£ = /3 S + ia s (i.e., a s is of the 
smallest value within all poles), we find 



Po 



(52) 



where a p = qo<x s > 0. This scaling form is verified in 
Fig. [2] (a): All the data from different systems char- 
acterized by distinct elastic constants (i.e., different B x 
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FIG. 2. (Color online) (a) Pinning strength po as a function 
of interface thickness £ or reduced temperature e (inset), for 
different values of elastic modulus B x . The result indicates a 
universal relation given in Eq. (|52[) : po ~ exp(— a p £), with 
Qf p = 0.6620±0.0008. (b) Change of interface thickness £ with 
the reduced temperature e. Note that larger value of e corre- 
sponds to lower system temperature; also for comparison, the 
numerical value of atomic layer spacing is ao = 2n/qo ~ 6.28. 



values) can be scaled onto a single universal curve obey- 
ing Eq. (|52|) (except for very small values of po < 10~ 14 
for which numerical errors would be too large), where 
a p = 0.6620 ± 0.0008 as determined from data fitting. 

Note that numerical results in Fig. [2] (a) seem to im- 
ply an asymptotic behavior of po ~~ ^ as the system ap- 
proaches the melting point (i.e., e — > 0). However, since 
this is a subcritical bifurcation system (with hexagonal 
symmetry), the interface thickness £ remains finite and 
the pinning strength po would never vanish as e — > (30| 
in both real systems and the full PFC model. 

A universal scaling behavior can be also identified for 
the surface tension er, although the form of scaling is 
different. As shown in Fig. [3] where the results are calcu- 
lated from Eq. (HTj) . values of a for systems of different 



FIG. 3. (Color online) Surface tension a as a function of £ or e 
(inset). All data for different values of B x collapses on a single 
universal curve of a vs. £, which can be fitted into 2 power- 
law relations: a ~ £-2.62±o.o2 £ or sma jj interface width and 
a ~ £-2.969±o.oo5 £ Qr more diffuse interfaces. The crossover 
occurs at an intermediate width value £ ~ 33. 



elastic modulus B x are well fitted into a single scaling 
curve of a vs. £. This data collapse works well for all 
range of interface thickness £ in our calculations, yield- 
ing a power law behavior a ~ £- a <* ; although with two 
power-law exponents found in two distinct regimes: For 
thin enough interface a a = 2.62 ± 0.02, while more dif- 
fuse interface results in a faster decay of a determined 
by OL a = 2.969 ± 0.005. The crossover between these two 
scaling regimes is identified in Fig. [3] Note that gener- 
ally surface tension a becomes larger for larger value of e 
(lower temperature), and also for sharper interface with 
smaller B x (see the inset of Fig. [3]) , as expected. 

Another important parameter governing interface mo- 
tion is the kinetic coefficient Co, which determines the 
relation between the interface velocity and the thermo- 
dynamic dri ving force and is of large interest in solidifica- 
tion studies [3J, [35| . Note that the expression of Co given 
in Eq. (|46|) . if neglecting the last average density term, 
is similar to that determined by Mikheev and Chernov 
from classical density functional theory [34[ which can 
well describe recent results of molecular dynamics simu- 
lations [35[ . The anisotropic feature of kinetic coefficient 
identified in previous studies has also been incorporated 
in Eq. (|46|). as the amplitude profiles (A® and Uq) vary 
with the orientation of liquid-solid interface. 

On the other hand, here we focus on a system different 
from these previous studies [although the general form 
of interface equations (|44|) . (|49"]) , and (|48[) is applicable 
to both cases] : Instead of using interface undercooling as 
the thermodynamic driving force [3J, l35| , here we study 
the isothermal solidification process in pure materials, 
and the driving force originates from the supersatura- 
tion in atomic density at uniform temperature. We can 
then examine the kinetic coefficient for different isother- 
mal systems, each with a specific e value and the cor- 
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FIG. 4. (Color online) Kinetic coefficient £b as a function of 
reduced temperature e, for different values of B x . An estimate 
of Co ~ e 2 is shown for comparison. Note that the thermo- 
dynamic driving force here is not interface undercooling, but 
the density supersaturation in isothermal systems. 



FIG. 5. (Color online) Propagation of interface position h 
with time t, as evaluated from analytic solutions (|57[) for 
\F \ = l.lpo and ((59j for \F \ = 0.9p , with e = 0.3 and 
B x = 10. Positions of h are shown in unit of atomic layer 
spacing ao = 2n/qo. 



responding liquid-solid coexistence conditions. Results 
evaluated from Eq. (|46)) are given in Fig. |4l showing the 
increase of Co with decreasing temperature (i.e., larger e). 
Data of different elastic modulus B x fall around a power- 
law relation Co ~ e 2 , although with large deviation found 
at small B x and e values (i.e., B x — 10 and e < 0.02). 



B. Planar interface dynamics and pinning 



If neglecting the lattice pinning effect, the dynam- 
ics of interface profile is trivial: h(t) — h(Q) + vot, 
with a constant interface growth rate vq = -Fb/Co = 
— <7qA?1o(^o + A)/Co- However, as will be shown below 
the lattice coupling effect plays a significant role in the 
description of interface dynamics, even for the simplest 
case of layer-by-layer growth considered here. 

The dynamical equation governing a planar interface 
profile h(t) is derived from Eq. ([50]). i.e., 



In the case of planar growth along the normal y di- 
rection, the lateral variations in the interface equations 
of motion can be neglected, resulting in an effective ID 
system. In the frame co-moving with the interface at a 
velocity vq, the outer equation (|49|l in the liquid region 
(y > 0, and y — » y — Dot) leads to a steady-state form 



v d y Sn + (qq + A) 



dySn 



0, 



(53) 



where /3o = — e + 3uq (+oo) — 2gno(+co). The liquid- 
state chemical potential variation is given by 8jJt(y > 
0) = (Qo + M8no, satisfying a far-field boundary con- 
dition Sn(+oo) = A which represents an external growth 
condition of constant flux coming from the liquid bound- 
ary. In the solid side (y < 0), the outer solution yields a 
constant S/j,(y < 0) = <5/i(0). From the interface condi- 
tion vqAuq — —d y 5[i\ + as determined by Eq. (|35)> and 
the continuity of 5/i, we obtain the steady-state solution 



g f<5/zoexp(- 
[ + A, 



y< o, 



y >o 



where = (q^ + /3o)A?t,q. Thus from Eq. 
effective driving force is given by 

F Q = A = -g 2 An°(^ + A). 



(54) 



3 the 



(55) 



CoVn = (odh/dt = F - p sin(q h + f)+r] v , (56) 

which can be solved exactly in the absence of the noise 
term r] v , as given in the following. 



1. \Fo\ > po: continuous growth mode 

When the magnitude of external driving force exceeds 
the lattice pinning strength po, the exact solution of the 
interface position is written as 



h(t) =± — arctan < (3 tan 
qo { 



go 
2Co 



{F 2 



p 2 o) 1/2 t±r 



(57) 



with "+" for F > and "-" for F < 0, where (3 = 
[(F -po)/(F +p )} 1 / 2 and r is determined by initial 
condition, i.e., 



f 1 

t = arctan < — tan 



*M0) + £ - - 

2 W 2 4. 



(58) 



This situation could occur at small enough e (i.e., high 
temperature growth) and diffuse enough interface, and 
hence small enough pinning force (see Fig. [5]), given a 
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certain flux condition Fq . Despite its form of continuum 
description, the solution (|57[) yields the jumps of distance 
2tt/<7o (= ao) for interface position, which is exactly one 
discrete lattice spacing along the y direction of growth. 
As shown in Fig. [5] which gives the numerical evaluation 
of Eq. (|57| , the liquid-solid interface propagates "contin- 
uously" due to the overcoming of lattice pinning, while 
the discrete lattice effect can still be preserved in this 
continuous growth mode as a form of growing steps of cto 
spacing. In this case the average velocity of interface can 
be calculated as v = (dh/dt) = ±(F 2 - pl) 1/2 /( - 



2. \Fq\ <po-' activated/nucleated growth mode 

For the growth condition of lower temperature and 
sharper interface (with larger pinning strength) or weaker 
driving force such that \Fo\ < po, the exact solution of 
Eq. (|55|) without noise is different: 



h(t) = — arctan 



1 + 



l-^/o(«) 



tp- 



/*>, (59) 



where ft' = [(po - 
F$y/H/Co}, and 



F )/(po + Fo)} 1 / 2 , f (t) = exp[q (p 2 Q 



_ tan [f h(0) + $-$\-£ 
T ° tan[f/i(0) + f -f] +P'' 



The interface growth rate is then given by 



dh/dt 



(4/Co)(Po-foK/oW 
[1- W)f+/3'2 [l + r / (t)f 



(61) 



At large time t ^> 1, v n = dh/dt — > 0; the interface is thus 
locked/pinned by the underlying crystalline potential at 
a position h — —(2 arctan j3' + <p — ir/2)/qo (satisfying 
sm(qoh + (p) = Fq/pq). This pinning phenomenon is il- 
lustrated in Fig. [5] which shows the numerical evaluation 
of the analytic solution ([591 . 

Thermal fluctuations should then play an important 
role on the process of lattice growth and interface moving, 
and the full stochastic dynamic equation (1561) with noise 
term r\ v governed by Eq. (|4"T)) should be used. This 
will become a stochastic, escape problem in a potential 
system [30| . and the liquid-solid front would propagate 
via an activated process to overcome the pinned lattice 
site, a procedure analogous to thermal nucleation. To 
illustrate this depinning process, we rewrite Eq. (|56[) as 



dh/dt 



dUh 
' dh 



Vh, 



(62) 



where the effective potential Uh = — [(po/<Zo) cos(qoh + 
tp) + F h}/^ Q and the noise rjh satisfies = 
2D S(t-t'), with D Q = D/d = tfq 2 r Q k B T/( . From the 
corresponding Fokkcr-Planck equation we can determine 
the Kramers' escape rate 



t 2tt 



d 2 u h 



dh 2 



d 2 U h 



dh 2 



,-AU h /Do 



(63) 



which represents the rate of an atom hopping/escaping 
from a metastable lattice site "a" (determined as a lo- 
cal minimum of potential Uh) to a nearest lattice site 
with lower potential, via overcoming a potential barrier 
AUh = Uh(b) — Uh(a) where "6" indicates the top loca- 
tion of the barrier (i.e., a local maximum of Uh). In Eq. 
(1531 r is the escape time of atoms. It can be shown that 



d 2 U h /dh 2 



-d 2 U h /dh 2 \ b = q (p 2 



Fg^/Co, and 



the potential barrier 



AC4 



<7oCo 



(Pi 



F 2 o 



1/2 



|.F | arccos(|-F |/po) 



(64) 

(which is always positive for \Fq\ < Po)- 

In this thermally activated process, the lattice nucle- 
ation growth rate is given by I = aoR = ao/r, where 
do = 2ir/qo is the spacing of atomic layers along the 
growth direction y as shown in Fig. [T]and Sec. HU From 
Eqs. (|6"3"|) , ([6^1) and the expression of D , we obtain the 
standard Arrhenius form for thermal nucleation: 



where 



I = I e- E ^ kBT , 



h = {pl-F 2 f/ 2 /^ 



(65) 



(66) 



(60) and the activation energy E a is determined by 



E a 



(pI f 2 



1/2 



|F |arccos(|F |/po) 



It is important to note that in the general form of Eq. 
(|65l). both the prefactor I and activation energy E a are 
actually dependent on temperature (e) and also elastic 
constants (B x ), as can be seen from their expressions in 
Eqs. ([51]) and ||6~T|) . For a simple example, if setting a 
growth condition of Fo = a^po (|«o| < 1) f° r all temper- 
atures or e values, we have the rescaled activation energy 
E' a - E a m = (2/gg)[(l - a 2 ) 1 / 2 - \a \ arccos |a |]po, 
showing the same behavior of temperature and interface 
width dependence as that of po (see Fig. [2|) . On the other 
hand, considering constant driving force Fo at different 
temperatures would lead to more complicated tempera- 
ture and width dependence of E a and Iq, as shown in 
Fig. [5] The results there are obtained from numerical 
evaluations of Eqs. ||6"T|) and (|6"6"]) . Fig. [5] (a) shows 
that both E' a and Jo increase with e (i.e., the decrease 
of temperature). At low temperatures with large e and 
small interface thickness £, p 3> F and Eq. (|57|) yields 
E' a cx po ■ We would then expect the activation energy to 
follow a universal scaling relation similar to that of pin- 
ning strength po: E' a ~ exp(— a p £). A deviation would 
occur for large enough ^ (i.e., small enough e and high 
enough temperature) due to similar order of magnitudes 
between p and F values, as has been verified in Fig. [S] 
(b). Also interestingly, all the numerical data of prefac- 
tor Iq for different values of elastic modulus (B x ) is found 
to collapse on a universal curve Iq ~ exp(— where 
aj = 0.54 ± 0.01 as obtained from data fitting. All these 
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FIG. 6. (Color online) The rescaled activation energy E' a = 
EafiTo and the prefactor Jo of the nucleation growth rate as a 
function of (a) reduced temperature e and (b) interface width 
£, for different values of B 1 and a constant Fo = 10 . In (b), 
all data of Io are fitted into a scaling form Io ~ exp(— a/£), 
where cti = 0.54±0.01; for the data of E' a , the scaling relation 
of po [~ exp(— a p £)] is shown for comparison. 



results indicate that in the PFC model the temperature 
dependence of nucleation rate / is not exactly Arrhenius, 
but shows a more complicated nonlinear behavior. 



V. DISCUSSION AND CONCLUSIONS 

We have constructed a nonadiabatic amplitude rep- 
resentation and examined the sharp-interface limit of 
the 2D single-component PFC model. Our main find- 
ings include the coupling and interaction between meso- 
scopic description (for "slow" variation of structural am- 
plitudes) and microscopic scales (for "fast" variation of 
the underlying crystalline structure), and also the re- 
sulting interface pinning effects which are incorporated 
in a generalized Gibbs-Thomson relation for interface 
dynamics. The strength of the corresponding pinning 



force, and also the value of solid surface tension, have 
been found to obey universal scaling relations with re- 
spect to the liquid-solid interface width. Temperature 
dependence of the interface kinetic coefficient has also 
been examined for isothermal systems driven by density 
supersaturation. The scale coupling effects have been il- 
lustrated in an example of planar interface growth, which 
shows the crossover between two distinct growth modes, a 
"continuous" or nonactivated mode for high-temperature 
or strongly-driven growth and a "nucleated" mode for 
low/moderate temperature or weakly-driven growth, as 
a result of the competition between the external ther- 
modynamic driving force and the lattice pinning/locking 
effect. Such thermal nucleation process in the growth 
regime of \Fq\ < po is analogous to a pinning-depinning 
transition in the absence of quenched disorder. 

Note that the sample application given in Sec. IIV Bl 
can be viewed as the lowest order approximation of inter- 
face growth, and the scenario of nonactivated/continuous 
vs. activated growth modes identified is qualitatively 
consistent with the classical crystal growth theory of 
Cahn [36[ based on the concept of critical driving force, 
and also with that found in early work of roughening 
transition by Chui and Weeks @ and Nozieres and Gal- 
let although more complicated parameter depen- 
dence on temperature and material elastic property is 
determined here. More general study should involve non- 
planar surface/interface evolution and dynamics, so that 
details of dynamic roughening and faceting transition in- 
volving this new parameter dependence can be obtained. 
Also, all the above derivations can be readily extended 
to (2+l)D PFC growth systems with 3D bcc or fee lat- 
tice structure. The corresponding interface equations are 
expected to be of the same form as Eqs. l[15|). P5|) . and 
the generalized Gibbs-Thomson relation (J3U); the latter 
could be also reduced to a driven sine-Gordon form (i.e., 
the form of the (2+l)D Hwa-Kardar-Paczuski equation 
fllj ] ) governing the interface profile h(r,t): 



(odh/dt = F + aV 2 h + - \Vh\ 2 - p Q sm{q h + tp) 



■ Vv, 

although with more complicated expression of the co- 
efficients. It is important to note that generally Fq 
and A are functions of h and lateral coordinate r in 
both (1+1)D and (2+l)D cases, as can be seen from 
S/i(0, s) = Sfj,(y = h,x) in Eq. (|35]), Thus the above 
driven sine-Gordon form is actually a much more com- 
plex nonlinear equation of h, and is closely coupled to the 
conservation condition P5)l and the outer equation P^)l 
which reflect the conserved dynamics of atomic number 
density. On the other hand, if neglecting the spatial de- 
pendence of Fq and A at lowest order approximation [i.e., 
approximating Fq and A by the planar constant result Eq. 
(|55l) for the case of weak surface fluctuations] , the Hwa- 
Kardar-Paczuski equation with spatially-constant coeffi- 
cients can be recovered. In further studies it would be 
interesting to identify the properties of the corresponding 
dynamic roughening transition as compared to the pre- 
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vious renormalization-group results [T(| [H| (noting that 
from our derivation coefficients in Eq. (|68p are generally 
temperature dependent), although one would expect such 
transition to be near the lowest order result of |Fo| = po- 

Another important topic to be addressed, based on 
the nonadiabatic amplitude representation and sharp- 
interface approach developed here, is the anisotropy 
along different crystal growth directions, particularly for 
surface tension (<j), kinetic coefficient (Co)) an d pinning 
strength (po). This would yield detailed properties of 
facet formation and anisotropic, orientation-dependent 
roughening transition. Note that all the results given in 
Sec. Illll and llVl are for interfaces oriented along the y di- 
rection. For other orientations we expect the same forms 
of interface equations of motion due to similar deriva- 
tion procedure, but with different values/expressions of 
coefficients. For example, as shown explicitly in Eqs. 
(TTl"]) - ([i"4")) and Fig. [U atomic layers oriented along y or 
(±V / 3i + y)/2 direction are not equivalent to those along 
x or (±x+V3y)/2 direction, due to different layer spacing 



(a vs. a /^/3) and hence different degree of scale cou- 
pling and pinning effect. Similar properties are expected 
for 3D PFC models of various symmetries, although with 
more complicated results anticipated. 
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Appendix A: Derivation of lattice pinning term in 
the interface equation of motion 

To derive the interface equation from the nonconserved 
dynamic equations for Aj, we perform integration of 

J_j*dUdjjAj* on Eq. (|2"5)) and also summation over all 
the resulting equations for Aj and A* . Corresponding to 
the last scale- coupling term in Eq. (|28|). we get 



j 

a iq h 



, \ r u + x y 

(duAf) / dU'f Pj (A°, n°) e^u'+H) + , 



A 



U+Xy 



U+Xy 



dA° / dU'f; x e llioU ' + / dAf / dU'f P2 e^ oU ' + dA° 3 dU' f* p3 e^ oU ' 



U+Xy 



+ c.c. 



a iq a h 



dUA^ u [P pi {U + X y ) - f* pi (U)] + / dUA° 2 *e^ u [f P2 (U + X y ) - f P2 (U)] 



+ / dUA° 3 e^ u [f* P3 (U + X y ) - f; 3 (U)} + 



c.c. 



iq ° h / dU e ii!oU 



Alduf^+ATduf^+Aldu^ 



Pa 



(Al) 



where f Pj (U) = f Pj (A®(U), n[j(C7)), and we have used 
q X y = 4tt and [f p AU + X y ) - f P ,(U)]/X v ~ duf } 



iPj- 



It 



is then straightforward to show that (|A1|) is equivalent 
to the lattice pinning term po sin(<7o^ + f) appearing in 
Eq. (|3ip . with the pinning strength po and phase ip de- 
termined by Eq. (|33|) . 



Appendix B: Matching between inner and outer 
regions and the Gibbs surface condition 



In the inner region, the solution of Eq. AMI) for Green's 
functions G satisfying the corresponding boundary con- 



ditions at U, U' = 0, ±C has been given in Ref. [If, i.e., 

r+lTT <!■ TP Q'\ - J - U ' 5 ^ S - 0<U' <U<( 

[,,,) ~\-U5(S-S'), 0<U<U'<( 

C-(TT TT> ?>\ - I US ( S - 5 ')< ~t<U' <U<Q 
u \u,o,u ,o ) - i u ig(g_ S i^ -( <U <U' <0. 

(Bl) 

Substituting solution (|B1|) into Eq. ([35]) leads to 



v(S)e F dU'[n° - n° (+C)] = h{U, S) - Ai(0+, S) 
Jo+ 

-U (duhk + £ / du '[^ ~ Vo(+0l (B2) 

-v(S)e / dU'[n° - ng(-C)] - fii(U, S) - £i((T,S) 
Jo- 

-U (%Mi)_c + £ / dU 'tio - Vo(-Q], (B3) 
Jo- 
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where Eq. (|B2|) applies to 0+ < U < ( and (|B3|) applies 

to-c<f/<o-. 

From the rescaling (U,S) — (it/£, es/£) in the inner 
region and (U°,S°) = (eu/£,es/£) for the outer region, 
the matching conditions (|43l) can be rewritten as 



(duh) u= ±£ = £ (%°£i u %°=± £ e ■ ( B4 ) 



Also, due to eC = e (/£ "C 1 we can carry out the expan- 
sion for the outer solution 

p,T*(0±,S°) = Ar(±eC, S°) T e( (du'fiT*)^ + 0{e 2 ). 

(B5) 

Evaluating Eq. dB2j) with U = +C and Eq. (|B3]) with 
U = — C, and using conditions (|B4[) and (|B5|) . we find 



- «£ 2 / C rff/'K - ng(+C)] = Ai ut (0+, S°) - Ai(0+, 5) 
Jo 

H [ C dU'[rjS - %"(+0] + 0(e 2 ), (B6) 
Jo 

ve J__dU'[n° - ng(-C)] = Ai ut (tr , S°) - £i((>- S) 
-£ /^K " %"(-C)] + 0(e 2 )- (B7) 



Adding Eqs. (|B6|) and (|B7p. using the Gibbs surface 
condition (|42[). and considering £ = £/£ ^> 1, we get 



Mi(0, 5) = Mi ut (0, 5°) + vC 2 / dC/K - n° (+oo)} + fj m , 

Jo 

(B8) 

where fj m = r) m /e = {J °° du[?7 I -^(+oo)]-/° oo du[fj%- 
f/ 1 (-oo)]}/2 (with ^ = ^/e). Rewriting Eq. (|B8j in 
the original scale (u, s) and substituting into the inter- 
face equation (|4l7|) , we can obtain the generalized Gibbs- 
Thomson relation given in Eq. ([44]) . 

To derive the standard interface continuity condi- 
tion (|4"8|) , we apply the matching condition (|B4|) to 
Eq. (|39p and expand the outer result around u = 0; 
that is, (dt/Ai)±c = ^(^Mr t )±ec = e(%"Mi ut )o± ± 
ff 2 C(9^ o /i? ut ) ± + C(e 3 ). Keeping terms up to 0(e) and 
neglecting the noise effect in the outer solution would 
then yield Eq. (|48| in the original scale. 

Note that the Gibbs surface condition can actually be 
determined from Eqs. (JB6J) and JB7]) up to (9(e). Sub- 
tracting Eq. (|B6[) from (|B7[) and neglecting the noise 
terms, given the continuity of /ii and jl° ut at u = we 
obtain 



C 



_dC/ [ng-ng(±0] ^C(e 2 ). 



(B9) 



C 



Returning to the original scale and noting £ ^> 1, at 0(e) 
we can recover Eq. (|42|) for the Gibbs surface. 
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